Formfactor perturbation expansions and 
confinement in the Ising field theory 



S. B. Rutkevich 

Institute of Solid State and Semiconductor Physics, SSPA "Scientific-Practical 
Materials Research Centre, NAS of Belarus", P. Brovka St. 17, 220072 Minsk, 
Belarus 



Abstract. We study the particle spectrum M n (h) in the two-dimensional 
ferromagnetic Ising field theory in a weak external magnetic field h. According to 
Wu and McCoy scenario of the weak confinement, pairs of fermions (domain walls) 
are coupled into bosonic kink-antikink bound states at small h > 0. Fluctuations 
with more than two fermions also contribute to the wave functions of the compound 
particles, leading to the multi-fermion corrections to their masses M n (h) in higher 
orders in h. We describe a perturbative procedure, which allows to account both 
multi-fermion fluctuations, and the long-range confining interaction between fermions, 
and leads to the formfactor expansions for the renormalized parameters of the model. 
We obtain integral representations for the third-order multi-fermion correction to the 
mass M n (h), which arise from the regular correction to the kernel of the Bethe-Salpeter 
equation. 
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1. Introduction 

In recent years much progress has been achieved in the understanding of the scaling 
limit of the two-dimensional Ising model, which is known as the Ising Field Theory 
(IFT), for a review see pQ. Providing direct information about the Ising universality 
class in two dimensions, IFT can be viewed also as a continuous dynamical model of the 
one-dimensional uniaxial ferromagnet. Being, perhaps, the simplest relativistic model 
describing confinement of topological excitations, IFT can give a deep insight into some 
nontrivial aspects of confinement in the particle and condense matter physics. 

IFT contains parameters m and h, which are proportional to the deviations of 
temperature T and magnetic field H from their critical values in the two-dimensional 
lattice Ising model, m ~ (T c — T), h ~ H . At the critical point m = 0, h = 0, IFT 
reduces to the Conformal Field Theory with central charge c = |, which Euclidean 
action Aqft describes free massless Majorana fermions. It has two relevant operators, 
the energy density e(x), and the order spin operator a(x). IFT can be defined as the 
perturbation of Ising Conformal Field Theory by these two operators, which is described 
by the action [2] 



In fact, only one dimensionless parameter rj = m/\h\ 8 ^ 15 determines the physics of IFT. 

IFT being not integrable for generic h and m, admits exact solutions along the 
directions h = and m — 0. The line h = 0, m ^ corresponds to Onsager's exact 
solution J3|. Fermions remain free here, but gain the mass \m\. In the disordered 
(paramagnetic) phase m < these fermions are ordinary particles, while in the ordered 
(ferromagnetic) phase m > they are interpreted as topological excitations (kinks), 
which separate regions with oppositely directed spontaneous magnetization. Nonzero 
magnetic field h > induces interaction between fermions, breaking integrability of IFT 
at m 7^ 0. On the other hand, IFT has a remarkable exact solution at m = 0, h ^ 
containing eight massive particles, which was found by A. B. Zamolodchikov [4]. 

Beyond the integrable directions, IFT can be studied by approximate methods - 
numerical and analytical. An effective numerical method known as Truncated Conformal 
Space Approach was invented by Yurov and Alexei B. Zamolodchikov [5], [6]. Fonseca 
and A. B. Zamolodchikov |7J modified this technique and applied it to analysis of 
analytical properties of the IFT free energy continued to complex values of the scaling 
parameter r\. 

For analytical study of IFT for h and m close to the integrable directions, it is 
natural to exploit perturbation expansions. Form-factor perturbation theory developed 
by Delfino, Mussardo and Simonetti [8] has been applied [8] , [9] to calculate the variation 
of the particle mass spectrum and the decay widths of non-stable particle for small 7], 
i.e. near the line m = 0. One could expect, that the perturbation expansion at m / 
and small h should be more simple, since (unperturbed) IFT is free at h = 0. Though 
this is really the case in the high-temperature phase m < 0, the small- h expansion at 




(1) 
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m > turns out to be rather non-trivial due to the long-range attractive potential 
between neighbouring fermions, which is induced by the external magnetic field h > 0. 
This attractive interaction can not be accounted by the straightforward formfactor 
perturbative theory at small values of h, and leads to confinement of fermions. 

The effect of a small magnetic field h, which brakes the Z^-symmetry in the ordered 
phase m > in IFT, can be qualitatively understood by the following simple arguments 
first developed by McCoy and Wu [10] . At h — 0, two ferromagnetic ground states 
| 0+) and | 0_) with spontaneous magnetizations +<r and —a have the same energy. 
A weak magnetic field h > removes degeneration decreasing the energy of the state 
| 0+), and increasing the energy of the state | 0_), which becomes metastable. In order 
to generate a domain of the metastable phase in the stable surrounding, one needs to 
add the energy proportional to the length of the domain. In other words, two domain 
walls bounding such a domain attract one another with the energy 2h a I proportional 
to their separation I, see figured! The long-range attraction leads to confinement: all 
domain walls are coupled into pairs at arbitrary small h > 0. Elementary excitations 
now are the domains bounded by two kinks, while an isolated kink gains infinite energy. 
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Figure 1. Two kinks interact with the energy 2hal. 



The mechanism of confinement outlined above is quite general in one-dimensional 
systems. It is realized in such continuous one-dimensional models as multi-frequency 
Sin-Gordon model [TTJ, 5- state Potts field theory [12], and in the discrete Ising spin 
chain [13]. Confinement of topological excitations in one-dimensional antiferromagnet 
has been observed experimentally by Kenzelmann et al. [14|. On the other hand, there 
is a lot of similarities between confinement in the IFT and in 't Hooft's model for two- 
dimensional multicolor QCD [15], see the discussion in [2j. Accordingly, the fermions 
and their bound states in IFT in the confinement regime are used to call as "quarks" 
and "mesons", respectively. 

At small h the weak confinement regime is realized in IFT. In this regime, the mass 
spectrum M n of mesons is dense in the segment [2m, 00). Two asymptotic expansions 
describe M n at h — > in different regions of this segment. Near the edge point 
2m (i.e. for fixed n at h — > 0) one can use the low energy expansion in fractional 
powers of the magnetic field [7], [2]. On the other hand, for n ^> 1 and h — > 0, the 
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semiclassical expansion in integer powers of h can be applied [2], [TB]. Derivation of the 
both expansions are based on the perturbative analysis of the Bethe-Salpeter equation, 
which determines the meson mass and wave function in the two-quark approximation. 
The latter implies, that one approximates the meson wave-function (the eigenstate of 
the IFT Hamiltonian) by the two-quark state, neglecting multi-quark (four-quark, six- 
quark, . . . ) contributions to it. The two-quark approximation is asymptotically exact 
in the limit h — > giving correct meson masses in the leading order in h. However, 
starting from the second order in h, it is necessary to take into account the virtual multi- 
quark fluctuations. Note, that multi-quark effects are essential also for such interesting 
phenomena as decay of unstable mesons, and inelastic meson scattering. 

The second order multi-quark correction to the meson mass was obtained by Fonseca 
and A.B. Zamolodchikov [2j. These authors demonstrated also, that the multi-quark 
corrections could come up in the weak-coupling expansions of the meson masses M n in 
three ways: 

(i) through the radiative corrections of the quark mass and self energy; 

(ii) by renormalization of the long-range attractive force between the neighbouring 
quarks (the "string-tension"); 

(iii) by modifying the regular part of the Bethe-Salpeter kernel, which is responsible for 
the pair interaction between quarks at short distances. 

It turns out, that only the first contribution (E]) gives rise to the second-order correction 
to the meson mass, while (Jn]) and (TTIT|) should show up only in the third order correction, 
which is still unknown. 

Extension of the weak-coupling expansions for the meson masses to the third order 
in the magnetic field presents an interesting problem, which we address in this work. It 
could give us some insight into the role of the multi-particle fluctuations in the composite 
particles in non-integrable models exhibiting confinement. Since multi-quark effect are 
responsible also for the decay of unstable mesons, this should manifest itself in some 
form in the perturbative meson mass spectrum near and above the stability threshold. 
Note, that an accurate numerical calculation of the lowest meson masses was reported 
in reference [2], which clearly indicates contribution of the multi-quark fluctuations. 

Since the problem outlined above is rather involved, here we shall concentrate 
only on three parts of it. First, we extend the semiclassical expansion of the original 
(written in the two-quark approximation) Bethe-Salpeter equation to the third order 
in h. Second, describe the formfactor perturbative technique, which is suitable to deal 
with the multi-particle fluctuations in systems with confinement. Finally, we obtain the 
integral representations for the "local" multi-quark correction of the meson masses, i.e. 
corrections (lull) induced by renormalization of the local interaction between quarks. 

The paper is organized as follows. In section [2] we describe definition of IFT 
and its operator content. In two subsequent sections we summarize briefly the recent 
progress in the theory of the weak confinement in IFT: section |3] introduces the Bethe- 
Salpeter equation and its weak-coupling expansions, and section 0] contains preliminary 
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discussion of the multi-quark corrections to the meson masses. In section we develop 
a formfactor perturbative procedure, which is modified to a system with a long-range 
confining interaction between fermions. It is based on the partial diagonalization of the 
Hamiltonian in the fermionic number, and allows one to effectively account the multi- 
quark fluctuations by " dressing" the fermionic operators. In section [6] we describe a 
compact integral representation for the local third order correction to the meson mass, 



which is analysed further in |Appendix B[ |Appendix A| contains perturbative solution 



of the "bare" Bethe-Salpeter equation to the third order in h. Concluding remarks are 
presented in section [7J 

2. The model 

Ising field theory is the Euclidean field theory, which describes the scaling limit of the 
two-dimensional lattice Ising model in the critical region T — > T c , H — > 0. It is defined 
by the action 

1 Z* 00 r 1 /*00 

Atft = 7T~ / ipdip + ipdip + im^ d 2 x — h / d 2 x a(x). (2) 

^ J -oo L * J -oo 

Here x denotes a point in the plane M? with cartesian coordinates (x(x),y(x)) and the 
complex coordinate z = x + iy, d = ^(d x — id y ), d = ^(d x + id y ). Action (j2J) is covariant 
under rotation, and becomes Lorentz covariant after the Wick turn y — > it. 
Corresponding to action (J2j) Hamiltonian can be written in the form 

H = H + hV } (3) 
where H = — u(p) a f (p) a(p), 

J-oo 27r 



V=- dxa(x 



anda;(p) = (p 2 +m 2 ) 1 ^ 2 is the spectrum of free fermions. Fermionic operators (p'), a(p) 
obey the canonical anticommutational relations 

{a(p),at(p')} = 2ir6{p-p'), {a(p),a(p')} = {a f (p), at(p')} = 0. 

Commonly used are also fermionic operators a(/5), a\(3), corresponding to the rapidity 
variable (3 = arcsinh(p/m): 

a(P) = ^(p) 1/2 a(p), at(/3) = uip) 1 ' 2 a t(p). 

Notations 

- - - ,Pat> = a f (pi) . . .a+O^lO), (pi,...,pj\r| = (0|a(pi)...a(pjv), 
= a t (/3 1 )...a t (/? JV |0), = (0|a(/?i)...a(#v) 

for the fermionic basis states with definite momenta will be used. 

The order spin operator er(x) = cr(x, y)| y=0 in the ordered phase T < Tq (i.e. 
at m > 0) can be determined in the infinite line x e R as the normally ordered 
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exponent [TTlfTg] : 

a(x) = a : e^ x)/2 : , (4) 
P(x) 



/>oo 

/ dx'( X (x',y)«9 y x(x',y))| 3 
.7 x 



2 | \ 5 i/ / y/v \ 1 J / J ly=0 ' 

X(x,y) =i / ^^=fat(-p)e^)^-a(p)e-^^ 



2vr ^/^y 

where a = m 

l/8 2 l/12 e -l/8 A 3/2 ig the 

zero-field vacuum expectation value of the order 
field (spontaneous magnetization), A = 1.28243... is Glaisher's constant. 

Alternatively, operators <r(x) can be completely characterized by their formfactors 
(Pi, . . . , PkW(0)\P[, . . . ,P' N ), which explicit expressions are well known [T9], [7J. In the 
ordered phase 

(p 1 ,...,p K \a(0)\p' 1 ,...,p' N ) (5) 
= i (K +N)/ 2- -Q tanh(^) J] tanh(^Z^) ]J coth f^-S\ , 

0<i<j<K ^ ' 0<k<q<N ^ 0< <K 

0<t<N 

if (K + N) is even, and (ft, . . . , (3 K \a(0)\(3[ } . . . , /3' N ) = for odd (if + iV). The right- 
hand side in §5§ contains factors coth y lf* J , which are singular at at (3 S = (3' t . These 
kinematic singularities should be understood in the sense of the Cauchy principal value 



coth 



2/4 



Note, that the Wick expansion holds for formfactors ([51) of the spin operator. For 
example, 

a(ft,/5 2 |a(0)|ft,^) = (Pi\o-(0)\(3' 2 )(f] 2 \o-(Q)\(3[) 
-(PiW(0M)(P 2 \a(0M) + </? 1 ,/3 2 |cx(0)|0>(0| t r(0)|^,^>. 

3. Bethe-Salpeter equation 

The meson energy spectra AE n (P) can be formally determined from the eigenvalue 
problem: 

H | $ n (P)) = [AE n (P) + £ vac ] | $ n (P)>, (6) 
H = Ho + hV, 

P\$ n {P)) = P\$ n {P)), (7) 
where P is the total momentum operator, 



A r°° dp t . . . . 

J-oo 27r 



-oo 

and E VSbC is the ground state energy, which is proportional to the length of the system 
L. 
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The eigenvalue problem ([2D is quite difficult, since the Hamiltonian contains the 
order spin operator cr(x), which is highly nonlinear in fermionic fields. A significant 
simplification can be provided by the two-quark approximation [7j, [20]. It implies that 
one replaces the exact Hamiltonian eigenvalue problem ((ED, © by its projection to the 
two-quark subspace F 2 of the Fock space F: 

V 2 H | £ n (p)) = [AE n (P) + E vac ] I $ n (P)>, (8) 
P|$ n (P)) = P|$ n (P)>, 
I $„(P)> e F 2 . 

Here V n denotes the orthogonal projector onto the n-quark subspace ¥ n of F. Tildes 
distinguish solutions of ((§} from those of the exact eigenvalue problem (jBJ). 
In the momentum representation, equation (jSJ) takes the form [2] 



cu(P/2 + p)+^(P/2-p)-AP(P) *p(p) = / 



HA- 

G P (p\k)* P (k)^ } 
2n 



where j denotes the Cauchy principal value integral, 

(P'/2 - p, P'/2 + p|$(P)) = 2u5{P' - P) * P (p), 

G P (p|A;) = £(P/2 + p, P/2 - p|P/2 + k, P/2 - fc), 

1/4 



G(pi,p 2 \k u k 2 ) = —{P2,pi \a 



h, k 2 ) 



u(p l )uj(p 2 )uj(k l )uj(k 2 )\ 1 /' 2 
u(pi) + uj{k 2 ) oj{p 2 ) + u(ki) w(pi) + oj[ki) uj(p 2 ) + uj{k 2 ) 



(9) 

(10) 
(11) 



Pi 



+ 



P2 - h 

Pi ~P2 



Pi 

ki - k 2 



P2 



Uj(pi) + Uj{p 2 ) u(k\) + U>(k 2 ) 

and /o = 2/icx = A m 2 is the "bare string tension" . Index n is omitted in ([9]) 
that typ{p) is an odd function of p, and 



M- Note, 



Gp(p|*) 



(p - k) 2 (p + fc) 2 



where Gp (p|fe) is regular at real p and g. The pole terms in Gp(p|/c) produce after 
the Fourier transform the long-range linear attractive potential /ol x l proportional to the 
distance |x| between the two quarks. The regular term G^p C9 \p\k) is responsible for the 
local interaction between quarks vanishing at the distances ^> vrC 1 . 

Equation (Q is the Bethe-Salpeter equation written in a generic momentum frame. 
It simplifies in two cases. 

• In the frame of the centre of mass of two quarks [7j, P = pi + p 2 = 0: 



00 dk 

2co(p)-AE(0) * (p) = /of TT" 



%{k) 



2tt 2u{p)u(k) 



(12) 



uo(p) + uo(k) 
p — k 



+ 



p 



2u;(p)cj(£;) 
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In the infinite momentum frame (see Appendix A in [2]), P 
^ 2 M 2 



m 



u- 



*(u) = fo 



F(u\v) $(v 



i 



2tt 



where the scaled variables u = (p± — P2)/ P and v 
and 



(ft 



F(u|u) = [(1 - w 2 )(l - 
$(u) = lim * P (Pu), 



-1/2 



1 



OO: 

(13) 

q2)/P have been used, 
1 + Wf 



UV 

+ T 



The following large-P asymptotic behaviour of AE(P) was assumed in 
deriving (fl3"j) from (jUJ): 

M 2 



AP(P) = |P| 



2 P 



0(|P 



-3> 



111 



(14) 



Bethe-Salpeter equation ([9]) and its particular cases f[T2"j) . f[T3"j) are the linear singular 
integral equations |21j. Different techniques [2], [7], [16] have been been developed 
for their perturbative solutions in the weak-coupling limit A — > 0. Fonseca and 
A.B. Zamolodchikov calculated [2] several initial terms in the low-energy expansion 
(for fixed n and A — > 0) for eigenvalues of equation (fi~3~]) 



M 2 

n 

4m 2 



2, 



34 57_ 

5 V 175 + 280 
1894^ 23983 2 2 



t 10 + 



23 £ w 
7875 
3313 z 



1543 z n 
12600 



t 8 + 



13 



1120 7T 



(15) 



3031875 242550 J 10080 tt 
where t = A 1//3 , and (— z n ) is the zero of the Airy function, Ai(— z n ) = 0. The leading 
term in the above expansion reproduces the old result of McCoy and Wu [10] . 

To the second order in A, semiclassical expansions (for n ^> 1 and A — » 0) for M 2 , 
and for AE(0) were found in references [2], and [16], respectively. We extend the former 
expansion to the third order in A using the technique, which was applied previously in 
the similar discrete-chain problem [13]. This calculation is described in Appendix A 
the result reads as 



M n ,2a 

— \ = cosh 2 e n , 
where 6 n solves equation 

sinh 29 n - 26 n = 

and 

Si{e n ) 1 / 5 
s 2 (9 n ) 



2n\{n - 1/4) + 2\ 2 S 1 {6 n ) + 2\ 3 S 2 {6 n ) + 0(A 4 ) 
sinh 2 9 n ' 



1 

12 



1 

4 cosh 2 6 n 



6 



sinh(2# n ) V 24 sinh 2 9 n 

' -{-999 n - 3fl n [648cosh(2fl n ) + 228cosh(4# n ) 



192vr sinh 6 {26) 

+56 cosh(60„) + 15 cosh(8# n )] + 546 sinh(2fl n ) + 363 sinh(40„ 
+ 170sinh(6fl n ) + 33sinh(8#„) + sinh(12# n )}. 
To the second order, (fT6l) - (lT9l agrees with [2J. 



(16) 
(17) 

(18) 
(19) 
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4. Beyond the two-quark approximation 

Eigenvalues AE n (P) of the Bethe-Salpeter equation are not the same as the 
eigenvalues AE n (P) of the initial problem ([6]): 

AE n (P) = AE n (P) + 5E n (P). (20) 

The difference SE n (P) is caused by the multi-quark corrections, which are ignored in 
(Q, but contribute to AE n (P). The exact meson energy spectra should have the form 

AE n (P) = {Ml + P 2 ) 1/2 . (21) 

due to the Lorentz covariance of IFT, but this form does not hold [2 J for the meson 
energies AE n (P) determined in the two-quark approximation. 

In the P -> oo limit equation ([20]) yields due to ([14]) and ( 12~T1) : 

M 2 = M 2 + 5M 2 , 

where 

8Ml = lim [2P5E n (P)\. 

P—*oo 

The first analysis of the multi-quark corrections to the meson masses has been 
done by Fonseca and A.B. Zamolodchikov [2]. They claim, that multi-quark corrections 
treated perturbatively in A should modify the Bethe-Salpeter equation ([9]) to the form 

/°° dk 
G P (p\k)* P (k)—. (22) 

Here e(p) and / are the renormalized quark dispersion law and the renormalized string 
tension, respectively. The renormalized kernel Gp(p\k) is assumed to have the structure 

G P (p\k) = G P (p\k) + AG { ; eg) (p\k), (23) 

where Gp(p\k) is the original kernel ([TT]) . and the correction term AG p" e9 \p\k) = 0(A), 
being regular at k = ±p, effectively modifies the pair interaction between quarks at 
short distances < m _1 . 

Note, that the renormalized quark energy does not have the Lorentz covariant form 

e(jp) = (p 2 + m 2 ) 1/2 + 5e(p) = (p 2 + m 2 ) 1/2 + Ae(p), 

since quarks are not free particles at h > due to their confinement. Assuming 
Ae(p) = 0(|p|~ 3 ) at p — > oo, one can define the "dressed" quark mass m q from the 
large-p asymptotics of e(p): 

2 

771 

e(p) = \p\ + ^ + 0(\pn. 

There are no nonperturbative definitions of renormalized quantities in equation ( [221 . 
Instead, it is expected, that they can be determined order by order by their power series 
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in A: 

m 2 q = m 2 (1 + a 2 \ 2 + a 3 A 3 + •••), (24a) 

5e(p) = 5 2 e(p) + 5 3 s(p) + 0(A 4 ), (246) 

/ = / (l + c 2 A 2 + c 4 A 4 + ---) (24c) 

AGp 9) (p\k) = A 1 G i ; eg \p\k) + A 2 Gp 9 \p\k) + 0(A 3 ). (24d) 

Let us summarize briefly, what is known about the coefficients in the above 
expansions. Fonseca an A.B. Zamolodchikov [22] analyzed the exact integral 
representation for the coefficient a 2 in (124 a\) , and obtained from it the value 

a 2 = 0.071010809... (25) 

On the other hand, one can expand a 2 into the sum 

of the second order (in A) diagrams with three, five, . . . , quarks in the intermediate 
state. Contribution of three-quark diagrams into a 2 was estimated in reference [7] 

a 2>3 ~ 0.07... (26) 

We obtain its exact value 

° 2 - s = li + l2^= - 07094 "" (27) 
this calculation will be presented elsewhere. Comparison of fl27|) with fl25l) shows, that 
the second order radiative correction to the quark mass is essentially determined by the 
three-quark contribution. Diagrams with five and more quarks in the intermediate state 
give less than 0.1 % of a 2 . 

The term of order A 2 in expansion (12461) for 5e(p) was found by Fonseca and 
A.B. Zamolodchikov [2]: 

S,e(p) = - X -^-. (28) 

yF! 2 w(p) 8 u 5 {p) K ' 

They have given also strong arguments, that coefficients c 2 k in expansion (124 cp should 
be simply related with coefficients gj in the well known weak-h expansion [19] for the 
vacuum energy E vac , 

E vac = Lm 2 (- l - A + g 2 A 2 + g 3 A 3 + g 4 A 4 + . . . J , (29) 



namely 



c 2 fc = -2g 2k+x . (30) 



In particular, c 2 = —0.003889 . . . 

It is not difficult to modify the weak coupling expansions (both low-energy and 
semiclassical) to account renormalized quantities in the Bethe-Salpeter equation fl22|) . 
and to the express multi-quark correction 5M n in terms of coefficients in (I24ap -( l24~dj) . 
It turns out [2], that for calculation of the meson masses M n to the third order in A, 
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it would be sufficient to know the renormalized quark mass m q and the string tension 
/ to the third order in A, and the "regular" term AGp (p\q) in fl23l) to the linear 
order in A in the limit P — > oo. To this end, one needs to determine two unknown 
quantities: the third order correction to the quark mass (coefficient 03 in (124 al) ). and 
the kernel AiG^ 9 \p\k) in (124 dl) . In fact, we need only the diagonal part of the latter, 
A 1 Gt 9 \ P \p). 

The problem of explicit calculation of a 3 and A 1 G P " e9 ' ) (p\k) is quite difficult. Here we 
do not try to find its complete solution. Instead, in subsequent sections we shall obtain 
several representations for these quantities in terms of formfactors of spin operators o~(x) 
and their products a(xi)cr(x2). 



5. Diagonalization of the Hamiltonian in the fermionic number 

Bethe-Salpeter equation ([9]) is approximate, since the IFT Hamiltonian ([3]) does not 
conserve the number of fermions - the "bare" quarks. Let us try to find a unitary 
operator U(h), which transfers operators generating "bare" fermions into operators 
generating such "dressed" fermions, that their number would be conserved by the 
evolution operator. It is clear, that the two-fermion Bethe-Salpeter equation, written 
for these "dressed" fermions, should be exact, and it could be identified with the 
renormalized Bethe-Salpeter equation (1221) . 

Let at(p), a(p) be the set of creation/annihilation operators of the "dressed" 
fermions, which are related with the "bare" ones by the unitary transform 

a(p) = U{h) a(p) Ufa)' 1 , a f (p) = U{h) a f (p) U(h)- 1 

with operator U(h) depending on the magnetic field h. We shall also underline all 
"dressed" operators and states: 

A = U(h)' 1 AU(h), |$) = f/(/ l )- 1 |$). 

Expanding U(h) into the power series in h 

00 

U(h) = l + J2h n Fn, 

71=1 

we obtain the set of equalities following from the unitarity condition U(h) U(hy = 1: 
T\+T\ = 0, 

^ 2 = ^ + A, A f = -A, 

^3 = ^±^ + K, yt = _y, 

Denote by N_ the operator of the number of " dressed" fermions 
K= r^at(p)a(p), 

./-ex, 27r 
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and by V_ n the projector operators onto the subspaces of n "dressed" fermions. For an 
operator A acting in the Fock space let us separate the diagonal and off-diagonal parts 
in the "dressed" fermion number n, A = Ad + A s , where 

oo 

A d = ^2V n AV n , and A s = A - A d . 

n=0 

We require, that Ti and N_ commute, [H, N] = 0, or, equivalently, 

H s = 0. (31) 

Rewriting fl3Tj) as 

H s = (U{h)WJ{h)- l ) s = ^ 

one obtains 

(^l + hT 1 + h 2 T 2 + h i T- i + ...){H + hV){l + h^l + h 2 Fl + h 3 fl + ...)) =0. (32) 
Let us collect linear in h terms in (1321): 



(p\V\k) 

(j>\Fx\h)= ( (LA f ^ n{p)^n{k). (33) 

From here on we use compact notations \k) = \ki, Aw*;)); (p\ — {Pn{p)i ■■■iVi\i 
ui{p) = lj(pi) + ... + oj{p n {p))i an d so on - Equation f[3"3"j) defines (T\) S) but does not 
impose restrictions on (Fi)d- We fix the latter by the condition (Fijd = 0. 
In the second order in h one finds from (|32|) : 

(^m + Uv^ - Tm^x + [A,K ] + = 0. (34) 

This equation defines A s . We put = 0, and insert the intermediate state 
decomposition 

,>oo n (l) 



1 f°° 

i = £ m\ - 1 o> <o i + E ^ / i<w • • • • • • .2n( ff) i n 

9 „(,)=! "W" • / -« j=l 

into ( 1341) . providing 

i r <pizi?x?izifc) 



2tt 



(p|A|fc)= * V 



-w(p)][w(g) -w(fc)] 



■ _ uj(p)+u;(k) 



v (p|Zlg)(g|Z|fc) v (p|Zlg)(g|Z|fc) ] r , 

4- Kg)-w(p)] + V Mff) -"(*)] J 1 J 



9 

n(9)=n(fc) n(g)=n(p) 



for n(p) n(k). 

Note, that one can drop all underlining in the right-hand sides of equations (1331 and 
fl35|) . since = Similarly, we put Y d = 0, since equation ([32]) (in the 

third order in h) determines Y s only. 

In the rest of this section we shall consider, how the Hamiltonian TC acts in the 
subspaces with zero, one, and two renormalized fermions. 
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5.1. Vacuum sector 

In the vacuum sector, one obtains from ( |33l) . ( |34l) the standard Rayleigh-Schrodinger 
expansion ( |29l) for the IFT ground state energy: 

E vac = (0|W|0) = (0|f/(/i)(K + fcZM^lQ) = /i(0|V|0) + 8 2 E vac + 5 3 E vac + +0(/i 4 ), 
where 

fc^4'E« (36) 

n( 9 )#0 

5.1?. One-fermion sector 

In the one-fermion sector n(p) = n(k) = 1, and we find 

(p\n\k) = 2n5(p - k) oj( p ) + h( P \v\k) + 5 2 (p\n\k) + 5 3 (p\n\k) + o(h 4 ), (37) 

where 



1 

+ 



2 ^ L w (?)- w (p) -<*>(*). 

™(<3)#™(p) 



(38) 



5 3 (p\H\k) = +^ X>|V|g>(g| WWM*^ [1 - 5 n(q)Mp) ][l - 5 nW)Mp) ] (39) 



2 

g,g' 



+ 



[u(p) - iu(q)} [uj(p) - uj(q')] [u(k) - uj(q)] [u(k) - u{q>)\ 

^n(g),n(p)[l ~ ^n(q'),n(p)] [1 _ ^n(q),n(p)] 8 n {q'),n(p) 



u(q) - u(q') 



uj{q') — uj{p) uj(q) — ui{k) 



First, let us consider the linear term in h in the right-hand side of ( 1371) 

/■oo 

h(p\V\k) = -h dx(p|<7(x)|Jfe), (40) 



where 

i a exp[ix( k — p)} u(p) + u){k) 



{p\a(x)\k) 



p — k \uj{p)uj{k)] 1 / 2 
is the formfactor of the order spin operator (jl]) in the momentum basis. Integration in 
x in ( l4T)l) leads to the divergent result 

h(p\V\k) = -2m5(p - k)—^-r. (41) 

p — k 

This singularity is well known in the standard formfactor perturbative theory, where 
it appears as the divergency of the first order correction to the fermion mass, which is 
interpreted as a formal indication of confinement [8j [1] . 



Formfactor perturbation expansions and confinement in IFT 



14 



To give a meaning to equation (1411) . let us mention, that the generalized function 
5{q)/q is well defined and equivalent to —S'(q) in the class of the main functions 
ip(q) £ C 1 taking zero value at the origin, ip(0) = 0. So, one can formally write 

8(p - k)—^— = -5'(p -k) + C5{p - k) 
p — k 

with some indeterminate constant C. 

To get further insight, it is instructive to consider the matrix element (X.\hV\k), 
where the state (X| describes a "bare" quark located at the point X: 

<X| = ^ PX (P\- (42) 




Figure 2. One-fermion state represents a kink centred at X. 
For the matrix element of the order spin operator er(x) we get 

<xk(x)|*> = 



dp . _exp[rp(X — x) + ikx] oj(p) + io(h) 

— i cr : = 

2vr p-k ^uj(p)uj(k) 



sign(x-X)ae ifcX + iae ifcx 



dp exp[ip(X — x)] 



2tt 



P 



k 



x)] | 




1/4 


'iv(k)' 


1/4 | 








. W (P). 





(43) 



where $ denotes the Cauchy principal value integral. Here the first term in the right- 
hand side is non-local, while the second term is well localized near the diagonal x = X 
exponentially vanishing for |x — X| ^> mT x . Equation (1431) allows one to interpret the 
one-fermionic state (X| as a kink of width ~ m _1 centred at X, which devides the regions 
with magnetizations —o to the left, and +a to the right side of it, see figure [2j 
Substitution of ( 1431) into ( 140|) yields after integration in x 

/i(X|^|A;) = (/oX + C)e ifcX , 

where fo = 2hcf is the "bare string tension" , and the constant C is proportional to the 
length of the system L being infinite in the thermodynamic limit. Thus, Hamiltonian 
(J3J) acts in the one-particle subspace of "bare" quarks ¥ 1 as 

ViHVi = u(p) + f i + C, 

where x and p are the one-particle coordinate and momentum operators. The same 
formula written for "dressed" quarks 

V 1 HV 1 =e(p) + fx + C R (44) 
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gives us the perturbative definition of the renormalized quark dispersion law e(p) and 
renormalized string tension /. In the momentum representation (jUJ) takes the form 

(p\H\k) = 2n5(p - k) [e(p) + C R ] + 2mf5'(p - k), 

which should be compared with (1371) order by order in h. 

The second order term in (l3Tj) determines the leading correction to the quark energy 
5 2 £ (p) in expansion ( 124 fell 

h(p\H\k) = 2tt5( P - k) [5 2 £(p) + 5 2 E vac ], 

where S 2 E vac is given by (1361) . Explicitly, it can be described either by the formfactor 
expansion following from (|38p 

S 2 e{p) = S 2t3 e(p) + <*2, 5 e(p) + ■■■, (45) 
h 2 f°° dqi . . . dq n 5(q l + . . . + q n - p) 

•lim{p|cr(0)|gi, . . .,q n )(q n , . . . , gi|cr(0)|fc), 
fe— >p 

or by the equivalent integral representation 

/oo /*oo 
dx / dy|im(p| < 7(x,y)(l-7'iM0,0)|A:), (47) 
00 Jo k ~^P 

where cr(x, y) = exp(— ixP + jHo)a(0) exp(ixP — yTCo). 

Representations (I4"5"l) - (l4"7j) were first obtained and studied by Fonseca and 
A.B. Zamolodchikov [22j [2], we quoted their results in section 0] [see equations (125]) . 
(1261) . (1281) ]. We determine the exact large-p asymptotics of the integral (]4"6l) for n = 3 

<^(p) = ^ (A + 7^ V 0(P- 2 ), (48) 



2p V 16 12?r2 / 
which leads to (l2~Tj) . 

The third order term (|39l) in (1371) contributes both to the string tension /, and to 
the quark energy e(p). It determines S 3 e(p) and the constant a 3 in expansion (124 al) for 
the quark mass itlr. Calculation of a 3 would be of much interest for interpreting recent 
numerical calculation of masses of lightest mesons, see figure 7 in [2]. 

5.3. Two-fermion sector 

In the two-particle sector of "dressed" quarks, Hamiltonian acts as 

( P2,Pi \'H\ k 1 , k 2 ) = (2tt) 2 [u;(pi) + u(p 2 ))[d(p 1 - k 1 )S(p 2 - k 2 ) - 8{p 1 - k 2 )5(p 2 - h)) 

+h(p\V\k) + 5 2 (p\H\k) + 5 3 (p\H\k), +0(h A ) (49) 

where S 2 {p\H\k) and 63(p\H\k) are given by equations ( 1381) . ( 1391) with n(p) = n(k) = 2. 
Two initial terms in the right-hand side of (j4*9l) give rise to the "bare" Bethe-Salpeter 
equation ([9]). 

The explicit form of the second-order correction is 

S 2 (p 2 jh\H\ki 1 k 2 } = -4tt/ 5(pi +p 2 -fa- k 2 )5 2 Q(pi, p 2 \fa,k 2 ), 
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where 



<$20(pi,Pa|*i, ^ = ^^2^1^ d (^) 2j d -?^ (Pi +P2 - 9i - ■ ■ • - g 2i ) (50) 



3=2 

(p2,piW(o)\qi, • • • , i2j)(q2j, • • • , ?ik(o)|fci, ^2) 
1 

+ 



o;(gi) + . . . + u(q 2 j) - u(pi) - io{p 2 ) u{q x ) + . . . + uj(q 2j ) - uj(h) - uj(k 2 )_ 

Application of the Wick expansion to the formfactors in the integrand brakes (l50l) 
into the sum of diagrams. Some of them contain one or two products of the form 

{Pa or g 7 (g 7 <7 ) fcp) = ■ . P P- , 51 

MPaMqj)} 1 ' 2 [u}(kpp{q 7 )]V 2 p a -qp kp- qp 

which have two kinematic singularities in the integration variable g 7 at g 7 = p a , and 
g 7 = kp. Here V- denotes the "principal value" generalized function, 

11/1 1 
V- 



z 2 \z + iO z - i0 y 
Let us rewrite the singular factor in the right-hand side of (15T1) as 

' = P f„ - q L-a) + ^ 5{Va " W(Pq " ^' (52) 

Fa 1/7 1/7 iPa ( / 7 Jl"'/3 q~{ I 

where 

1 1 r 1 1 

V = 1 

(Pa - ? 7 ) ( fc /3 - <? 7 ) 2 [ (p a - ? 7 - iO) - Q7 - iO) (p a - <? 7 + iO) - g 7 + iO) 

Substitution of (1321 into factors (|5T|) leads to splitting of diagrams containing (one 
or two) such factors into several (two or four) terms. The resulting diagrams can be 
classified by the number of 5-functions 5(p a — kp) arising from the second term in the 
right-hand side of (1521 . 

(i) Diagrams with two (^-functions give rise to the vacuum energy correction 5 2 E vac . 

(ii) Diagrams with one 5-function contribute to the corrections S 2 e(pi) , 5 2 e (p 2 ) to the 
energies of two quarks. 

(iii) The rest diagrams are regular at p a = kp for a,f3 = 1,2. We denote their sum by 
8 2 Q^ re9 \pi, p 2 \k\, k 2 ). It determines (to the linear order in h) the kernel AGp (p\q) 
in the renormalized Bethe-Salpeter equation (122]) : 

A 1 Gp 9) {p\k) = § 2 g {re9 \P/2 + p, P/2 - p\P/2 + k, P/2 - k). 
6. Local multi-quark corrections to the meson masses 

It is not difficult to account perturbatively the regular correction term A 1 Gp £9 \p\k) in 
the Bethe-Salpeter equation ( 1221) both in the low-energy and semiclassical expansions. 
The resulting local multi-quark correction to the meson energy reads as 



S 3 E n (P) 



fju\P/2) & 
Am 2 dp 2 



AMp 9 \p\p) (53) 

p=0 
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in the low-energy case n ~ 1, and 



SA(P) = - 2 A lC g-> ( ^ az^) (54) 

|(pi-j? 2 )v| V 2 2 / 

in the semiclassical case n> 1. Here momenta p± and P2 are the solutions of equations 

Pl+P2 = P, 

ufa) + ufa) = (P 2 + Mlf'\ 

and 

_ Pi _ P2 
^(Pl) W(p2) ' 

To obtain the local multi-quark corrections to the meson masses, we rewrite ( |53l) . 
(|54|) in the rapidity variables (3\ = arcsinh(pi/m), f3 2 = arcsinh(p 2 /^), and then proceed 
to the limit P —>■ oo, yielding 



5 3 M 2 A 3 / 9 9 ^ 2 
— hm 



m 2 ^,^) 

m 2 8 (9^00 ^(9/?! <9/3 2 7 /3i=/3 2 ff 2 

in the low-energy case, and 

m 2 " M 2 -4m 2 j™ a 2 W 

in the semiclassical case. Here r] = arccosh[M n /(2m)], and 

W(p 1 ,p 2 ) = ^( Pl )u;(p 2 ) A 1 Gp 9 \p 1 ,p 2 ), (57) 

where pj = msinh/^- for j = 1,2. Function W(j3i,j3 2 ) determined by (J57jl admits a 
compact integral representation, analogous to (j!7j) : 

/>oo />oo 

H/(A,/3 2 )= / dx / dy hm (/? 2 ,/5i|a(x,y)(l-P -P 2 )a(0,0)|^,/? 2 ). (58) 



In Appendix B we extract from this function the leading at (/3i + /3 2 ) - * oo "irreducible" 
part W irr (Pi, (3 2 ) which determines 5 3 M%. It is expressed there in terms of the two- 
fermion matrix elements of the order spin operator pairs, which are explicitly known 
[22]. 

The third-order term 8^{p\H\^) in ( )49j) also contains regular and singular parts. The 
former contributes to the meson masses only in the forth order. The latter renormalizes 
the quark dispersion law and the string tension, which give rise to the " nonlocal" third- 
order multi-quark correction to the meson energy E n (P). It is expected [2], however, 
that in the limit P — > oo the nonlocal multi-quark corrections to E n (P) can be absorbed 
by renormalization of the meson mass and string tension. Thus, the third-order multi- 
quark correction to the meson masses should be described by representations (|55|) and 
( I56p . in which the "bare" parameters should be replaced by the renormalized ones 



m — > m q , X — > X R , where Xr = f /m q , and m q and / are given by expansions (124 aft . 
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7. Conclusions 

This paper is devoted to extension to the third order of the weak-h expansions for the 
meson masses M n (h) in the ferromagnetic IFT. There are four third-order contributions 
to it. The first one comes from the Bethe-Salpeter equation ( fl3l) written in the two- 
quark approximation. For the semiclassical expansion this contribution is described by 
equations (|T6l) - (fl9l) . for the low-energy expansion was already determined in reference 
[2], see (fl5l) . Three other contributions to M n (h) are due to the multi-fermion 
fluctuations. The local contribution arises from the regular radiative correction to the 
Bethe-Salpeter kernel. For this contribution, we obtain the integral representations ( 1551) - 
( |58|) . which are compact and appropriate for analytical analysis, and representation 
( IB. 161) . (1B.15j) . (IB. 141) . which we plan to use in the future numerical calculations. 



The two last multi-quark contributions to the meson mass come from the third order 
corrections 5^e{p) and 5^f to the quark self energy and string tension, which are 
contained implicitly in the formfactor expansion ( 1391) . Explicit extraction of S^elp) 
and S^f from (139|) is rather involved. Whereas for 5%f the result is essentially known 
[see (124 d) . (|29l) . (J30l)]. explicit calculation of the third order correction to the quark self 
energy and quark mass remains an open problem. 
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Appendix A. Perturbative solution of the Bethe-Salpeter equation 

In this Appendix we describe perturbative solution of the "bare" Bethe-Salpeter 
equation (fl3l) in the infinite momentum frame P — > 00 to the third order in the magnetic 
field h. 

Appendix A.l. Some exact relations 

It is convenient to rewrite equation (Tl3l) in new notations 

±1 \ x.i \n 2\-V2 2 M 2 -Am 2 2ha 

M 2 M 2 
Since 4>(—u) = — 4>(u), equation (fT3]) takes the form: 



1 c\d 

(u 2 - = p + -<f>(v) ^ + 4^-^ , (A.l) 



-1 



7T 



UV 1 — UV ' 

u 4 

2 {u-v) 2 



or equivalent ly 

(u 2 - u 2 U(u) = 2iapu + 4p \-u + (1 - u 2 ) d u ] £ ( A . 2 ) 

J_l 7T V — U 
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We shall require, that <j)(u) is a purely imaginary function in the interval (—1, 1). 
Set 

1 r 1 dv 



9{u) 



■4>{v) 



2,711 J_ 1 v — u 

for complex u ^ [—1,1]. Function g(u) is analytical in the region D shown in 



Im u 



lm u 



D 



Reu 



B 



Re u 



(a) 



(b) 



Figure Al. (a) D is the region of analyticity of g(u), (b) B is the region of analyticity 
of U(u), 



figure IXTT a). and 



2a 



g(u) = + 0(u ) at u — > oo, 

u 2 

providing 

2a = Res \ u=00 [g(u)u\. 

For real u G (—1, 1) we get 

00) = 9(u + iO) - g(u - iO), 
f dv (p{v) 

J _l 7T V — U 



(A.4) 



i[g(u + iO)+g(u-iO)}, 



and equation (1A.2j) takes the form: 

(u 2 -u 2 )[g(u + i0)-g(u-i0)] 

= 4ip [-u + (1 - u 2 ) d u ] [g{u + iO) + g{u - iO)] + 2\apu. 



(A.5) 



(A.6) 

Therefore, function 



Formf actor perturbation expansions and confinement in IFT 20 
Let us define two function in D: 

Ui(tt) = — 4ip[ — u + (1 — it 2 ) d u ]g(u) + {u 2 — u 2 )g{u) — iapu, 
U 2 {u) = Aip[-u + (1 - u 2 ) d u }g(u) + (u 2 - i> 2 )g(u) + iapu. 

Due to (IA.5I) . we have Ui(u + iO) = U 2 (m - iO) for u E (-1, 1 
U(iz) defined as 

Ui(m) for Imu > 0, 
U2(w) for lmu < 

can be analytically continued into the complex region B shown in figure lATV b). Note 
that U(iz) is even in B, and real in the interval (—1, 1). 

Let us solve differential equation flA.6j) with respect to g(u): 



U(u) 



9W 



l 

4~p 



dv 



U+(v) exp 



u — V 



[(1 -u 2 )(l -v 2 )} 1 / 2 



(l-u)(l + v) 
(1 + «)(!- 



n i(l-^ 2 )/(8p) 



(A.7) 



where XJ + (u) = Ui(tt) +iap-u, and the branch of the last factor in the integrand is fixed 
as 



arg 



1 - (u + iO) 
1 + (u + iO) 



arg 



1 + (v + iO) 
(1 - (v + i0)_ 



for real u, v G 



-1,1). 



Integration in (1A.7|) goes along the path C + shown in figure 

lm V 
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Function g(u) should be single- valued in D. The trivial monodromy behaviour of 
g(u) at u = oo is provided by flA.71) . if the following requirement is satisfied: 

,ni(l-^ 2 )/(8p) 



dw 



UJv) 



c 



(l-v 2 ) 1 / 2 



:i-v) 



cxp 



i v 



o. 



(A. 



where the integration path C is shown in figure IA2I The last condition determines the 
spectrum v n . If p -C 1, the integral in flA.81) is determined by contributions of the saddle 
points v = ±v of the function T(w ): 



T(w) 



1 - v 2 , /l + W ^ 

In l - w. 



(A.9) 



2 Vl-w, 

In the semiclassical limit n 3> 1, contributions of these two saddle points are well 
separated, and (1A.8I) yields the final asymptotical formula 



TK) / 1\ , / U+K + Aw) 

= pix \ n + P arg ( — ; 

4 V 4y» & \^l-(v n + Avy 



exp 



4p 



, (A.10) 



valid to all orders in p — > 0. Here 



AT(i/ + Aw) = T(z/ + Aw) - T(z/) - 



averaging (...) stands for 

/oo 
d Aw /(Aw) exp 
-00 

providing 

((Aw) 2 ^ +1 > = 0, 



iz/(Aw) 2 
4p(l - z/ 2 



z/(Aw) : 

1 - Z/2 



d Aw' exp 



ii/(Aw 



r\2 



4p(l - z/ 2 )_ 



(A.ll) 



(A.12) 



((Aw) 2 ^> = [4ip(l-z. 2 )/^r(l/2 + j)/r(l/2). 

At small p, Av, function U + (z/ + Aw) can be expanded as 

00 00 

u+(i/ + Aw) = 1 + J2 ca p l ( Au )' 

i=\ 1=0 

under appropriate normalization of <fi{u). 
Appendix A. 2. Perturbation Expansion 

To obtain the explicit semiclassical expansion for the spectrum v n , we need: 

(i) to calculate several initial terms in expansion flA. 12[) . 

(ii) to substitute (IA.12I) into flA.101) and to expand the expression in (...) in powers of 
Aw, 

(iii) to perform averaging in (lA.lOp by use of (lA.lip . 

(iv) to expand then arg(. . .) in flA.101) in powers of p. 
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Steps ©-([Iv]) are straightforward and well suitable for computer calculations, below we 
describe only the step (JI]). 

Let us write down the formal Neumann series solving equation ( 1A.1I) in the class of 
the generalized functions: 

<f>{u) = <p(°\u) + <pW(u) + 0(p 2 ), (A.13) 



0(°) ( u ) = l[S( u -v)-6(u + u)], 



1 



1 iu 2 — V 2 



(u 2 — V 2 ) 2 



+ 16- 



\-v 2 



u 2 — V 2 ) 3 



(A. 14) 



principal values are implied for the singular terms in ( 1A.14I) . Substitution of ( 1A.13I) into 
(JAU), HS3JI, USD yields: 

g(u) =^)( M )+^)( M ) + 0(p 2 ), 
a =a (o) +a (D +0 ( p 2 ); 

U(«) = U (0) (u) + U (1) (u) + 0(p 2 ), 



where 



9 {0) (u) 



9 9 ' 



U(0)( 



2tTZ/ 2 



2 -2-2v 2 + v\ (\-v 
2 + v 2 + In ' 



2v V1 + Z//J 

We skip lengthy expressions for g^ l \u) and U^(u). Note, that all singular terms 
at u = v and at u = —v cancel in U^(tt). The explicit expressions for the Taylor 
coefficients cu in (1A.12I) read as 

iu 

c io = - 7T> 



-6 + I0v 2 + 3v A - 3v e 



Cu 



Cl2 



< iixv' i {v 2 



2 + v\ n 

m 



-9 + 27u 2 - 37i/ 4 + 3i/ e v l - 3 



2vrz/ 4 

,6 ,,2 



In 



l + v 
1 - i/ 



l 

2' 



6ttz/ 4 (z/ 2 - l) 3 4vrz/ 5 *~Vl + i//' 

Imc2o = v cP^ I p. 

These constants are sufficient to obtain the equation determining v n to the third order 
in p: 



,5-6u 2 -9ut + 6^ 



4 ^ 1 4 W 12^(1 
r30 - 62u 2 + 54z/ 4 - 21^ + 18*/ - 3u, 



2 ) 

10 

n 



(A.15) 



+ 



10 - 4z/ 2 + 6z/ 4 + 4z/ 



" ^ln / 



47TZ/6 ~ V 1 + Z/ n 

In variables 9 n and A related with v n and p by 



+ 0{p A ). 



v n = tanh^n, 



A 



P 



4 cosh 9 r . 



equation flATBD takes form (jlSMigj). 
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Appendix B. Integral of the four-particle matrix element 

It was shown in section [6], that the local third-order multi-quark corrections to the meson 
mass can be expressed in terms of the integral of the four-particle matrix element 

/OO f-OO 
dx / dy lim (P' 2 ,P' 1 \a(x,y)(l-V -V 2 )a(0,0)\p 1 ,p 2 ) (B.l) 

over the Euclidean half-plane, see equations (|55|) . (1561) . In this section we extract from 
W(Pi,/3 2 ) the most important "irreducible" part W irr (Pi, p 2 ), and transform it to the 
form suitable for numerical calculations. 

It is straightforward to rewrite formula (15^1) . giving the local third-order 
semiclassical correction to the meson energy S^E^P), in terms of W(@i, j3 2 ): 

fha\ 3 4 m 6 Wljhjty , R9 , 

3 " ( } E n (P) (M* - 4m*) a 2 ■ (R2) 

Here 

P = m (sinh j3\ + sinh p 2 ), 
E n = m (cosh Pi + cosh p 2 ) 

are the meson momentum and energy, and Pi and p 2 are the rapidities of two quarks 
(forming the meson) at their collision, M n = (P 2 — P 2 ) 1 / 2 is the meson mass. 

The matrix element in the integrand in fIB.lj) can be expanded by use of the Wick 
rule 122]: 



G(P 2} Pi\PuP 2 )= lim <#,#|<7(x,y)<7(0,0)|A,&>= (B.3) 
G(Pi\Pi)G(p 2 \p 2 ) G(Pi\p 2 )G(p 2 \Pi) lfG(Pi,p 2 ) 



G G G\E(pi)E(p 2 ] 

Here we follow the notations of [22] 

x = (x, y) = (r cos 9, r sin 9), 

G(r) = (0\a(x)a(0)\0), 

G(r,9;Pi,p 2 ) = {0\o-(x)a(0)\Pi,p 2 ), 

(P'\a(x)a(0)\P) = 2tt5(P' - p) + G(r, 9; p'\p), 
imr 



E(r, 9; p) = exp 



sinh(/3 + it 



2 

and drop the explicit indication of position dependence for the matrix elements. The 
equality 

<A,AW«Mo)|o> CI'W'OJIA.A) 



[E(r,0;fa)E(r,e;fo)Y 

has been used in deriving (IB. 31) . Explicit expressions for the matrix elements G(r), 
G(r,9; Pi, p 2 ), G(r,9; P'\P) in terms of the solutions <p(r), x(r) of the Sinh-Gordon 
equation and associated Lax functions ^ r +(r, 9; P), \l/_(r, 9; 0) are known from |22j. 
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Proceeding to polar coordinates r, 9 in integral (IB.li . one can easily show, that 

/ d0G(r,0;ft,ft|ft,ft)= / d9G(r,9;(3 2 + (3,(3, + ftft + ftft + (3) (B.4) 
Jo Jo 

for arbitrary real (3. Proof of (1B.4P is based on the relations 

G(r, 9; (3 2 , (3, |ft, (3 2 ) = G(r, 0; (3 2 + i6, (3, + i#|ft + iO, (3 2 + iO), (B.5) 
G((3 2 + 17T, ft + i7r|A + ivr,ft + m) = G(ft, ft|ft, ft), (B.6) 

which follow from (1B.3|) and the similar properties of functions G(ft,ft), G((3i\(3 2 ), see 
[2J. Equality (1B.4P means, that the integral in its left-hand side is Lorentz invariant. 
Unfortunately, the integral 



oc 



rdr / dflG(r,#;ft,ft|ft,ft) 
Jo 

diverges at large r. It becomes convergent after substraction of the "reducible part", 
see flBljl : 

G(r, 9- (3 2 , /3i|/?i, (3 2 ) G(r, 0; /3 2 , ft ft) (B.7) 
- lim (/? 2 A|a(r,0)(Po + ^2)^(O,O)|ft,ft). 

However, the second term in the right-hand side here does not satisfy the monodromy 
property like (1B.6I) . This means, that the local multi-quark correction flB.2j) to the 
meson energy is not Lorentz invariant by itself. We hope that the Lorentz invariance 
form of 5^E n {P) will be restored in the third order in h after picking up all the terms 
contributing to it, as it happens [2] for the second order term S 2 E n (P). 

At the moment it is helpful to extract the "Lorentz invariant" term from the 
reducible part in (1B.7|) . Namely, we shall subdivide it as follows: 

lim ((3' 2 ,(3[\a(r,9)(Vo + V 2 )a(0,0)\(3 u p 2 } = AG(r, ^; ft, ft |ft, ft) + 5G(r, ^; ft, ft |ft, ft) 

where the first term satisfies the required monodromy property 

AG(r, 9- ft + iTT, ft + ivr|ft + ivr, ft + m) = AG(r, B; ft, ft |ft, ft), 

while the function SG(r, 9; ft, ft |ft, ft) does not satisfy such a property, but the integral 



rdr / d#5GM;ft,ft|ft,ft) 
Jo 

converges at finite ft, ft, and vanishes in the infinite momentum limit (ft + ft)/2 — > 
+oo. 

Note, that function AG(ft,ft|ft,ft) is analogous to function S((3\(3) defined by 
Equation (5.13) in Page 20 of reference [22], whereas 5G(ft, ft |ft, ft) is analogous to 
the zig-zag diagram (b) in figure 3 in Page 19. 

Let us obtain explicit expressions for AG(r, 9; ft, ftjft, ft). 
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1. Vacuum sector. 

(faPiHr, 0)7V(0, 0)1/3!,^) = (B.8) 
a 2 tanh 2 — — — exp{— irm[sinh(/3i + id) + sinh(/3 2 + i0)]} = 



2 a 2 tanh 2 ^L_A C os{rm[sinh(/? 1 + i9) + sinh(/? 2 + i0)]} - 

a 2 tanh 2 ^ 1 - ^ 2 exp{irm[sinh(/3i + \6) + sinh(/3 2 + i#)]}- 

The first and the second terms in the right-hand side of ( 1B.8I) should be assigned to 
AG(/3 2 ,/3i|/3i,/? 2 ) and to 5G(f3 2 , /3i |/5i, #2), respectively. 
2. Two-quark sector. 

(8' 2 ,(3[\a(r,e)V 2 a(O,O)\0 u 8 2 ) = J™ ^^<ft, #l*(r,0)|i», Vl )( Vl , V2 \a(0, 0)\fa, 8 2 ). 



This can be splitted into five diagrams: 

*\ y* 

) V [ 



+ ~~ ^X" + ) f (B.9) 



+ + 



Contribution of two former diagrams into AG(/3 2 , fli\8i, /3 2 ) are: 



4^(A|A)^(/? 2 |/? 2 ), (B.10) 



a 2 



•iwiA)] 2 . (B-ll) 



where 

Z" 00 d?7 

RiPilfh) — o" 2 exp[— imr(sinh/9i + sinh/3 2 )/2] / — expfimr sinh^] (B.12) 

J-00 2vr 

■ coth coth ^ 2^ ^ '^ y2 co ^ ^ ^ ^ si n [ mr (sinh /3i — sinh f3 2 )/2] 

f°° dr?' rf — 81 rf — 8 2 

+cr 2 exp[imr(sinh 8\ + sinh f3 2 ) /2] / — exp[imr sinh rf] tanh tanh , 

2?r 2 2 

5(/3 x |/3i) = lim i?(A|/? 2 ). 

Note, that < \m(3j < it for j = 1, 2 is supposed in ( IB. 121) . It is easy to verify, that: 

R(j3 1 + m\p 2 + m) = R(j3 1 \i3 2 ) 

The third diagram in (1B.9I) is proportional to the function f^(t), which determines 
the well-known large-distance asymptotics of the Ising correlation function [23]. Its 
contribution to AG(/3 2 , /3 2 ) is 

— > — 2(j 2 cos{rm[sinh(/5i + i#) + sinh(/3 2 + i6*)]}tanh 2 — - ^ 2 f^ 2 \mr), 

f i2) (t) = ~{[Kf(t) - Kim 2 - tK (t) K x {t) + \Klit)}, 
where K (t), K\(f) are the MacDonald's functions. 
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The forth and the fifth diagrams are equal to one another and can be written as: 

1 + exp [— imr[sinh(ft + i9) + sinh(ft + i9)}} 



r ^ --2, ft - ft 

= — icr tann 



■T 2 (r, ft + iff, ft + iO) + a 2 tanh 2 ^ 2 - ^ 



(B.13) 



+i a 2 tanh ^ - ^ 2 exp {-imr [sinh (ft + iff) + sinh(ft + i60]}A(?~, ft + id, ft + id) 
+ ia 2 tanh hzll r 2 ( r , ft + iO, ft + i9), 

where 

%(r, ft, ft) = B 2 (r, ft, ft) + U 2 {r, ft, ft) + V 2 (r, ft, ft) + A(r, ft, ft), 

and 

S 2 (r,ft,ft) = -i 



e imr(sinh??i+sinhr;2) 



OO 

oo 



U 2 {r,fa,fa) = - 
V 2 (r,ft,ft)=- 
AO", ft, ft) = -ie 



drji dm 
(2vr) 2 

^2 imr (sinh /3 t +sinh rj 2 ) qq+J^ 

2tt 2 



»7i - A ,i V2 ~ P. 
coth 



■ tanh 



m - p 2 , 1P1-V2 

tanh . 



im r (sinh j?i+sinh/3 2 ) co +h 

2vr 2 



Vi- Pi , ,V1~P2 
tann . 



imr (sinh /3i+sinh fc) 



d?7i d?? 2 ; 



imr (sinh rji +sinh '^2 ) 



+ u Vi-Pi , 1V2-P2, u Vi~V2 

■ tanh tann tanh . 

2 2 2 

Here we again suppose < Imft < 7r for j = 1, 2. 
Note, that 

7^(r,ft + i7r, ft + i7r) = exp [-imr (sinh fa + sinh ft)] T 2 (r,fa,fa). 

The two former terms in (IB. 131) contribute to AG(ft, Pi\Pi, fa), while all the rest terms 
in (1B.13P should be assigned to 5G(p 2 , ft|ft, ft). Thus, the irreducible part of the 
two-particle matrix element takes the form: 

G irr (r, 9- p 2 , fa\fa,fa) = G(r, 9; fa, fa\fa, fa) - AG(r, 6; fa, fa \fa, fa) (B.14) 



G(ft|ft)G(ft|ft) S{fa\fa)S{fa\fa) 



G 



a 2 



G{fa\fa)G{fa\fa) [R(Pi\p 2 )f 



G 



a 2 



i / G(p u p 2 ) y 

G\E(fa)E(fa) + G >&>™ 



where 

C 2 (ft,ft) = C 2 {r, 9; fa, fa) = 2a 2 tanh 2 

^U„„t.2 Pi ~ P2 cosj^fginh^ + i0 ) + s i n h(/5 2 + i0)]}[l - f^ 2 \mr)} 



+2 cT tanh 
—Ho 2 tanh 



2 

ft -ft 



1 + exp {—imr [sinh (ft + iff) + sinh(ft + i9)}} 



T 2 (r,fa+i9,p 2 + i9). 
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Integration of (1B.14j) in r and 9 gives the irreducible part of the factor W(Pi,j3 2 ): 

poo pn 

W irr (/3 1 ,/3 2 ) = rdr d0GW(r,0;&,/3i|#L,ft). (B.15) 
Jo Jo 

Integral in r here is convergent for small enough \/3\ — /3 2 \, and 

W irr (j3 1 ,p 2 )=W irr (j3 1 + p,fo + p) 

for arbitrary (3. 

On the other hand, the integral 



/*00 f'TT 

/ rdr / 

Jo Jo 



d66G{r, 9; fa, Pi\Pi, fa) 



converges and vanishes in the infinite momentum frame. 

So, the local multi-quark contribution to the third-order meson mass correction 
takes the form: 

2 fha\ 3 8m 6 WW(ft,fl») , * 

5sM « = M^4^ ^ ' (R16) 

with W irr (Pi, (3 2 ) given by flB.151) . Three other third-order contributions to come 
from the two-fermion Bethe-Salpeter equation in the infinite momentum frame ( fl3l) . 
and from the quark mass and string tension renormalizations. 
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